The Astrophysical Journal, xxx: xxx-xxx, 2008 XXX xx 

© 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. 



00 
O 
O 
(N 



oo 



43 
6 



(N 
> 

in 

m 

(N 
(N 

cn 
o 
oo 
o 



x 



ON THE DISTRIBUTION OF STELLAR MASSES IN GAMMA-RAY BURST HOST GALAXIES 1 

J. M. Castro Ceron 2 , M. J. Michalowski 2 , J. Hjorth 2 , D. Malesani 2 , 
J. Gorosabel 3 , D. Watson 2 , and J. P. U. Fynbo 2 

Received 2008 March 12; accepted 2008 XXX xx; published 2008 XXX xx 



ABSTRACT 

We analyse Spitzer images of 30 long-duration gamma-ray burst (GRB) host galaxies. We esti- 
mate their total stellar masses (Af*) based on the rest-frame iGband luminosities and constrain their 
unobscured star-formation rates (SFR) based on the rest-frame UV continuua. Further, we compute 
a mean M^/Ljc teB% = O.45M /L Q . We find that the hosts are low M + , star-forming systems. The 
median M* in our sample ((M*) = 10 9 ' 7 M Q ) is lower than that of "field" galaxies (e.g., Gemini Deep 
Deep Survey). The range spanned by M* is 10 7 M < M* < 10 11 M , while the range spanned by the 
unobscured SFR is 10 _2 MQyr _1 < SFR < lOM yr _1 . There is no evidence for intrinsic evolution 
in the distribution of M* with redshift. We show that extinction by dust must be present in at least 
25% of the GRB hosts in our sample, and suggest that this is a way to reconcile our lower, UV-based, 
specific SFR (<f> = SFR/M*) with previous claims that GRBs have some of the highest <f> values. We 
also examine the effect that the inability to resolve the star-forming regions in the hosts has on <f>. 
Subject headings: cosmology: observations — dust, extinction — galaxies: fundamental parameters 
— galaxies: ISM — gamma rays: bursts — infrared: galaxies 



1. INTRODUCTION 

It is central to contemporary cosmology to map the 
build-up of cosmic structure and star- formation (SF); 
and we know that the detection of a gamma-ray burst 
(GRB) is an indication that its host galaxy harbors mas- 
sive SF. GRBs are pulses of 7-rays from sources of cos- 
mological origins, and are the most luminous, photon- 
emitting events in the universe. As tracers of SF they 
have fundamental advantages: dust extinction has es- 
sentially no effect in their detection at 7-ray and X-ray 
wavelengths; GRBs can be observed to very high red- 
shifts; these redshifts can be measured from afterglow 
spectroscopy independently of the host magnitudes; and 
the selection of GRBs is independent of the host galaxy 
luminosity at any wavelength. That is to say that GRBs 
are unique eyes to gainfully look at the star-forming uni- 
verse. But the following critical questions should be an- 
swered in order to use GRBs as tracers of SF: What is 
the level of bias in optically-selected GRB host samples? 
And what is the intrinsic bias in the GRB-SF rate? 

A canonical model is well established for long-duration 
GRBs: they occ ur in star-fo r ming regions in star- 
forming galaxies (iBloom et a .1 120021: [G orosabe l et alJ 
l2003atlChristensen et al.ll2004 iFruchter et alJl2006h and 
are associated with stel lar core collapse events and hence 
with high-mass SF (e.g.[ Galama et al. 1998: Hiort fTet al.l 
20031 : iStanek et all 2003: Zch et al. 2004; CampanaZeFaD 



20061 ). The emerging picture, however, is complex. Most 



GRB host galaxies are faint and blue (|Fruchter et al.l 

1 This work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract 
with NASA. 

2 Dark Cosmology Centre, Niels Bohr Institute, University of 
Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Den- 
mark; josemari@dark-cosmology.dk, michal@dark-cosmology.dk, 
malesani@dark-cosmology.dk, darach@dark-cosmology.dk, 
jens@dark-cosmology.dk, jfynbo@dark-cosmology.dk. 

3 Instituto de Astrofi'sica de Andaluci'a (CSIC), c/. Camino 
Bajo de Huetor, 50, E-18.008 Granada, Spain; jgu@iaa.es. 



119991 : iLe Floc'h et al.ll2003r ). A few hosts show tenta- 
tive evidence of very high star-form ation rates (SFRs; 
ICharv et all 120021 : iBerger et alJl2003f ). but their optical 
properties do not appear typical of t he galaxies found in 
blind submillimeter galaxy surveys ijTanvir et al.l [20041 
IFruchter et aTll2006l ). 

It is currently debated how GRB hosts relate to other 
known populations of star-forming galaxies. At redshifts 
around 3 the UV luminosities of host galaxies and the 
metallicities of GRB sightlines are consistent with the 
expectation if hosts are drawn from the underlying popu- 
lation of all star-forming galaxie s weighted with the total 
SF density per lu minosity bin (Jak obsson et al.|[2 005 
iFvnbo et al.t l2008). With Spitzer 's (I Werner et all 1200 



IRAC (Infrared Array Camera; iFazio et all l2004n mid- 
infrared (MIR) photometry, together with optical and 
near-infrared (NIR) data, we can establish how the host 
galaxies relate to other star-forming populations in terms 
of total stellar mass (M*). This is essential if we are to 
understand the full range of properties of star-forming 
galaxies at high redshifts and fully exploit the potential 
o f GRBs as probes of c osmic SF. 

ICastro Ceron et all (|2006h studied a sample of 6 
long-duration GRB host galaxies observed with IRAC 
and MIPS (Mult iband Imager Photometer for Spitzer; 
iRieke et al.l I2004T) . They estimated their M* based on 
rest-frame K-band luminosity densities and constrained 
their SFRs based on the entire available spectral energy 
distribution (SED). In this work we extend the compu- 
tations to a sample of 30, but constrain only the unob- 
scured SFRs with the rest-frame UV continuum. This 
larger sample ought to allow for a more robust statisti- 
cal analysis, as well as to probe the distribution of M* 
in redshift space. To determine M+ we utilize rest-frame 
K flux densities (interpolated from observed IRAC and 
NIR fluxes). This extend s the data set presented by 
ICastro Ceron et ail (|2006h . yielding accurate values of 
M+ in a large host galaxy sample. To determine the 
unobscured SFRs we use rest-frame UV flux densities 
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(interpolated from observed optical fluxes). These unob- 
scured SFRs are lower limits to the total SFR of a galaxy 
due to the possible extinction by dust, and w e compare 
them with those of ICastro Ceron et al.l (|2006f ). Our pa- 
per is organized as follows: An overview of the sample 
selection is given in [J2] The analytic methodology is de- 
scribed in for the sample are derived in and 
$5] sees the computation of the unobscured SFRs. We 
conclude in SjHlwith analysis and discussion. We assume 
an Vt m = 0.3, fl\ — 0.7 cosmology with Hq — 70 km s _1 
Mpc -1 . 

2. DATA 

Our current sample is composed of 30 long-duration 
GRB host galaxies, t hree of them within the X-ray flash 
category (|Heisdl2003l ). We made the selection by requir- 
ing each host to have rest-frame if -band data available 
(for the purposes of this work we define K = 2.2 /mi 
± 0.3 /zm); thus the M ± estimator is well calibrated 
(|Glazebrook et ai1 l2004). An additional requirement for 
inclusion in the sample was the availability of the red- 
shift. 

This 30 GRB host galaxy sample spans a redshift in- 
terval < z < 2.7, with a median value z ~ 0.84. For 
comparison, the median redshift 4 of those GRBs de- 
tect ed prior to the start of operations of the Swift satel- 
lite (|Gehrels et al.ll2004h is (z) ~ 1.0, and that of those 
GRBs detected afterwards is (z) ~ 2.3; i.e., in this work 
we are chiefly looking at the lower end of the GRB red- 
shift distribution. Given the redshifts sampled, the rest- 
frame if-band data for 24 of the 30 host galaxies were 
obtained from the Spitzer Science Archive, where we ex- 
amined all publicly available hosts up to (and includ- 
ing) October 2007. The remaining 6 GRB hosts (980425, 
030329, 031203, 060218, 060505, and 060614) in the sam- 
ple have very low redshifts (z < 0.1), so in those cases 
^A oba ~ K\ leat (i-e., K\ ieat falls within the nominal width 
of K\ obs ) • The sample is presented in Table [TJ 

Each host (except GRBs 060218, 060505, and 060614) 
were imaged with IRAC. Detectors are 256 x 256 squared 
pixel arrays (scale = l'/2pixel _1 x l'/2 pixel -1 ; field of 
view = 5.21arcmin x 5.21arcmin). The instrumental 
PSFs (FWHM) are: channel 1 = T/66; channel 2 = l'/72; 
channel 3 = l'/88; channel 4 = l'/98. The optical and 
NIR data complementing IRAC in Table [1] were obtained 
from the literature. Two UV data points (GRBs 980425 
and 060505) come fr om our analysis of GALE X (Galaxy 
Evolution Explorer; iMartin et al"1l2003l I2005D data. 

3. METHODOLOGY 

For the MIR photometry we use official Spitzer 
Post Basic Calibrated D ata (Post-BCD) products (in 
ICastro Ceron et "all 120061 we carefully verified the Post- 
BCD with our own reductions). Host extraction is based 
on the archival imagery world coordinate system cali- 
bration and visually confirmed with optical and/or NIR 
comparison images from the literature. The median sep- 
aration between the host centroid in each IRAC image 
and the best set of coordinates published is well below 1". 
GRB 980425 is the only host galaxy resolved in the IRAC 
images a nd we have obtained it s photometry from the lit- 
erature (|Le Floc'h et al.ll2006D . None of the other GRB 
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log M. (M s ) [derived from K rest fluxes] 

Fig. 1. — Distribution of the total stellar mass (M 4 ) in GRB host 
galaxies. Filled histogram: 29 out of the 30 hosts in our sample, 
spanning a redshift interval < z < 1.5. We note that GRB 030429 
has been excluded from the histogram above and the calculation 
of the median M*. This is because its host was never detected 
at any wavelength and, consequently, no lower limit to M* can be 
estimated. The horizontal axis shows the inferred host M*, derived 
from interpolated rest-frame J\-band flux densities. The median 
M+ of the sample is (iW*) = 10 9 ' 7 Mq. For those host galaxies 
for which we have upper limits we estimate a conservative lower 
limit by extrapolating a flat spectrum (/„ oc v°) from the reddest 
NIR/optical detection (references in Table [T] column 4); then we 
split an area normalized to unit among the bins bracketed by the 
limits. For each GRB host for which we have detections we assume 
a normalized Gaussian distribution of the error bars in linear space. 
Then we allocate M* in proportion to the area of the Gaussian in 
each bin. Open histogram: Figure 2 from Savaglio et al. (2007) is 
shown for comparison. 



host galaxies of our sample are spatially resolved in the 
IRAC images, and their flux densities can be estimated 
using small circular aperture photometry. We measure 
the flux densities over a circled area of radius 2 pixels. 
In most cases this allows us to recover the emission of 
the host while avoiding contamination from other field 
sources located nearby. But in a few instances there was 
suspicion that the nearby field sources might be contam- 
inating our host galaxy photometry. As a sanity check 
we subtracted those field sources and redid the photom- 
etry. Field-source subtraction was performed using the 
detection output image given by the Source Extracto r 
software package (SExtractor: iBertin fc Arnoutsi ri996). 
where the detected sources were replaced by background 
noise. The photometry on the field-source-subtracted 
images was always consistent with the original aperture 
photometry. Aperture corrections have been applied to 
account for the extended size of the PSF. Based on the 
Spitzer Science Center recipe for Estimating Signal-To- 
Noise Ratio of a Point Source Measurement for IRAC 5 
we calculate conservative errors, including both statisti- 
cal and systematic estimates. Our flux density measure- 
ments and upper limits are given in Table [TJ column (3). 
We find that, of those hosts in our sample observed with 
channel 1, about 36% are detected. For channel 2 the 
rate is about 64%. T his is roughly of the sa me order as 
the detection rate bv lLe Floc'h et all (|2006h with IRAC 
channel 2 (44%), though we caution that both samples 
are incomplete and suffer from selection biases. 



1 http : //astro . ku . dk/~pall j a/GRBsample . html 



' http: //ssc . spitzer . caltech. edu/documents/irac_memo . txt 
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3 



TABLE 1 

Hosts: flux densities, total stellar masses, and star-formation rates 



Rcdshift IRAC K Icst (21 980 A) UV rcst (2 800 A) 



GRB Host 


2 


Ref. 


f,, in Tv ) 


Ch. 




Rcfs. 


n ^ ' 

|T0 9 Mr,) 


f,, ( li Tv) 


Rcfs. 


SFR(uv) 
(Mr vr~M 


(1) 


(2) 




(31 




(41 




V') 


(6) 




(7) 


970228 


0.70 


1 


<3.7 


1 


<4.2 


t, 30 


<5.7 


0.34 ± 0.16 


48 


0.60 ± 0.28 


970508 


0.83 


2 


<2.1 a 


2 


<1.8 


31, 30 


<3.5 


0.28 ± 0.15 


48 


0.71 ± 0.38 


970828 


0.96 


3 


3.9 ± 0.3 a 


2 


3.7 ± 0.3 


31, 3 


9.5 ± 0.9 


<0.44 


34 


<1.5 


980326 


~1.0 


4 


<2.7 


2 


<2.6 


t>30 


<7.1 


<0.015 


49 


<0.056 


980425 


0.0085 


5 


2 977 ± 101 


2 


6 389 ± 395 


t b , 32 


1.1 ± 0.1 


1 748 ± 173 


t c , 50 


0.39 ± 0.04 


980613 


1.10 


6 


38 ± l a 


2 


42 ± 1 


31, 6 


142 ± 3 


0.83 ± 0.11 


6, 30 


3.6 ± 0.5 


980703 


0.97 


7 


11 ± l a 


2 


11 ± 1 


31, 33 


29 ± 2 


3.2 ± 0.1 


~48 


10.9 ± 0.3 


981226 


1.11 


8 


4.5 ± 0.5 a 


2 


4.6 ± 0.5 


31, 8 


16 ± 2 


0.27 ± 0.03 


8 


1.2 ± 0.1 


990506 


1.31 


9 


2.0 ± 0.7 


2 


2.0 ± 0.8 


t, 34 


9.3 ± 3.8 


0.20 ± 0.04 


34 


1.2 ± 0.2 


990705 


0.84 


To 


19 ± l a 


2 


18 ± 1 


31, 10 


36 ± 2 


~1.8 ± 0.3 


10 


~4.7 ± 0.8 


000210 


0.85 


n 


3.3 ± 2.0 


2 


3.2 ± 1.8 


t, 35 


6.4 ± 3.6 


0.79 ± 0.07 


48 


2.1 ± 0.2 


000418 


1.12 


9 


4.8 ± 1.8 


2 


5.0 ± 1.9 


}. 36 


17 ± 7 


1.33 ± 0.04 


48 


6.1 ± 0.2 


000911 


1.06 


12 


<4.3 


2 


<4.3 


t, 37 


<13 


0.33 ± 0.08 


37 


1.4 ± 0.3 


010921 


0.45 


13 


11 ± 2 


1 


12 ± 2 


t,i3 


6.5 ± 0.9 


2.2 ± 0.1 


48 


1.6 ± 0.1 


020405 


0.69 


11 


<5.4 


1 


<5.3 


t, 38 


<7.0 


2.1 ± 0.1 


38 


3.7 ± 0.2 


020813 


1.26 


15 


<2.5 


2 


<2.6 


t, 38 


<11 


0.41 ± 0.08 


34, 38 


2.3 ± 0.5 


020819B 


0.41 


16 


97 ± 2 


1 


104 ± 7 


t.ifi 


47 ± 3 


4.3 ± 2.6 


16 


2.6 ± 1.5 


021211 


1.01 


17 


<2.2 


2 


<2.2 


t, 38 


<6.1 


0.20 ± 0.04 


38 


0.72 ± 0.15 


030328 


1.52 


18 


<29 


3 


<27 


t, 39 


<170 


0.56 ± 0.08 


39 


4.6 ± 0.6 


030329 


0.17 


19 


<4.9 


1 


<5.1 


t,40 d 


<0.37 


1.5 ± 0.2 


40 


0.14 ± 0.02 


030429 c 


2.66 


20 


<7.0 


4 


<7.3 


t, 20 


<124 


<0.060 


20 


<1.3 


030528 e 


0.78 


21 


<4.6 


1 


<3.8 


t.ii 


<6.5 


7.2 ± 1.4 


41 


16 ± 3 


031203 


0.11 


22 


216 ± 3 f 


1 


192 ± 13* 


t, 42 


5.3 ± 0.4 


119 ± 39 f 


51 


4.3 ± 1.4 


040924 


0.86 


23 


<2.9 


1 


<3.2 


t, 38 


<6.5 


<1.1 


38 


<2.9 


041006 


0.72 


24 


<2.9 


1 


<3.1 


}. 38 


<4.4 


<0.98 


38 


<1.8 


050223 


0.58 


25 


18 ± 2 


1 


18 ± 2 


t. 25 


17 ± 1 


<8.1 


25 


<10 


050525A 


0.61 


26 


<1.6 


1 


<1.6 


t, 43 


<1.6 


<0.48 


43 


<0.64 


060218 c 


0.03 


27 






20 ± 6 


44 


0.052 ± 0.015 


15 ± 3 


52 


0.053 ± 0.010 


060505 


0.09 g 


28 






298 ± 10 


45 


5.8 ± 0.2 


75 ± 6 


t c , 45 


1.9 ± 0.2 


060614 


0.13 


29 






3.8 ± 0.7 


46, 47 


0.15 ± 0.03 


0.37 ± 0.13 


53 


0.019 ± 0.006 



References. } Th is work: (fit) IBloom et alj|200ll Sit) IBloom et alj Il99l S3) IDiorgovski et al.l|200l| ; |H) IBloom et al.l [19991 ; 

S3) ITinnev et aLlll99Sl; (|6|) IDiorgovski et al.l 12003; i% IDiorgovski et a.1.1 119981; j8l) IChristensen et al.l 120051: J9T) [Bloom et al.l 12003; TToT) 
ILe Floc'h et al.112002}: lW) IPiro et alj|2002t lfl2l) [Price ; et al.|l2002al; (1131) IPrice et al.ll2002bl; 11141) IPrice i et alj|2003: 11151) IBarth et al.l|2003 ; 
[|16l)IJakobsson e t al. 2005bi n71) IVreeswiik et al.ll2003t (I18D iMaiorano et al.H2006t 11191) IHiorth et al. 2003; ( 20 )) IJakobflSon et alJI2004l; H21TI 
Rau et al. 2005; (22) Prochaska et al. 2004; (23) Wicrscma ct al. 2004; (24) Sodcrbcrg et al. 2006; (25) Pcllizza ct al. 2006; (26) Foley et al. 
120051; I27T i lPian et al.l 2006; (28) IColless et aUB OOl; (29) Delia Valle et alll2006a|; feojl llCharv et alj|200a <31i) ICagtro Ceron et al.112006 1; 
11321 ) Le Floc'h et al.ll200al33TTTVreeswiik et al.1119991; iM) ILe Floc'h et aJJI2003t 13&I IGorosabel et al.ll2003al; 113611 IGorosabel et al.M2003tj ; 
fl37nlMasetti et al.1120051: (1381) IWainwright et al.l ;2007; ( 39 ) IGorosabel et , al.H2005 at 11401) IGorosabel et al.H2005bl, (1411) IRau et al.H2004l ; (1421) 
IMalesani et al.112004: H43I1 1 [Delia Valle et al.ll2006bl; H4) IKocevski et al.H2007l ; {4Bh IThone et al 1 120081 ; (46) J. Hjorth 200 8, priv. comm.; 
ll47l)JCobb et alj|200d; 1481 1 IChristensen et al.1120041; l49l) IBloom et alj 120021" (50) Michalowski et al. 2008, in prep.; 1BH) IMargutti et al.l 
120071; (1521) ISollerman et alll200S 153ft IMangano et al.ll20OT 

Note. — Because host positions are well determined from previous broadband imaging, upper limits are quoted at the 2a level, 
while errors are lcr. All (UV, optical, NIR, and MIR) flux densities and magnitudes in this t able (includi ng those in the table notes) are 
corrected for foreground Galactic dust-extinctio n. Corrections to the IRAC wavebands follow Lutz (1999). For the UV, optical and NIR 
passbands we use the DIRBE/IRAS dust maps (Schlcgcl ct al. 1998). We adopt a Galactic dust-extinction curve A\/Ay, parameterized 
by R v = A V /E(B - V), with R v = 3.1 l|Cardelli et al.H19891 )' Col. (3); Our photometry of Spitzer's IRAC, publicly available, archival 
data O- Channel 1 = 3.6 jim; channel 2 = 4.5 /im; channel 3 = 5.8 fim; channel 4 = 8.0 fim. Col. (4): Interpolated flux densities for 
the rest-frame if -band ([[3J. The data used were obtained from these references. Col. (5): M* derived (Q from the rest-frame _R"-band 
flux densities listed in column (4), with M+/ LK ZBa t = OA Mq/ Lq. Col. (6): Interpolated flux densities for the rest-frame UV continuum 
(§3l. The data used were obtained from these references. Col. (7): Unobscured SFRs derived (J5)l from the rest-frame UV continuum flux 
densities listed in column (6). 

a Flux density values are taken from Castro Ccron ct al. (2006); we refine their error estimates. 

b Our photometry of 2MASS XSC Final Release (T wo Micron All Sky Survey Extended Source Catalog; released 25 March 2003; 
IJarrett et al . 2000; http : //www. ipac . caltech. edu/2mass/ ), NIR (Kg-ba,nd) archival data for galaxy ESO 184-G082 (/ 21 739 a = 6 510 fj.3y 
± 406/iJy). 

c Our photometry of GALEX (Galaxy Evolution Explorer; Martin ct al. 2003, 2005; http://galex.stsci.edu/ ), UV archival data for the 
host galaxies of GRBs 980425 (/ 2 267 A = 1 592 jtJy ± 162 jtJy) and 060505 (/ 2 2 &7 A = 72>Jy ± 10/xJy). 

d if-band is the closest passband, blueward of IRAC, for which this host has data available in the literature. It is a poorly constrained 
upper limit. We make use of it nevertheless, for methodological consistency (see §3). But we note that in this particular case, given the low 
rcdshift of the host, a much closer representation of reality is provid ed by the low e r limi t M* = 6.4 X 10 7 Mq (extrapolated from J-band 
and H-band data). This value is fully consistent with those cited by Thonc ct al. (2007) and references therein. 
e X-ray flash. 

f Because of the low G alactic latitu de (b = —4.6°) of this host we correct for dust-extinction overestimates. Following the recommendation 
bv lDutra et~aLl (2003) we scale the Schlcgcl ct al. (1998) reddening value multiplying by 0.75 and adopt Em w(B — V) = 0.78 mag. The 

UV flux density error (column 6) contains the additional 25% uncertainty estimated by Margutti ct al. (2007). 

s Redshift of the 2dFGRS Public Database (Two Degree Field Galaxy Redshift Survey; http://magnum.anu.edu.au/~TDFgg/), archival 
data for galaxy TGS173Z112. 
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Fig. 2. — Total stellar mass (M+) as a function of redshift for 
our 30 GRB host galaxy sample. Squares: Detections. Vertical 
bars: Distribution of the non-detection measurements of M*. The 
lower limits of these bars were calculated by extrapolating a flat 
spectrum (/„ oc u°) from the NIR/optical data (references in Table 
[T] column 4). 



For each GRB host we compute the flux density at the 
central wavelength of the rest-frame if-band by means 
of a linear interpolation in log space. We interpolate be- 
tween the IRAC channel and the closest passband, blue- 
ward of IRAC, for which data are available in the litera- 
ture. In a few cases the IRAC waveband corresponds to 
a rest-frame wavelength shorter than JiT-band, thus we 
extrapolate. The rest-frame X-band flux densities are 
shown in Table [JJ column (4), along with the appropri- 
ate references. In the same fashion we compu te, for each 
host, the rest- frame UV continuum (2 800 A; iKennicutti 
1998) flux density. We either interpolate between the 
two closest passbands, with data available from the lit- 
erature, that bracket the rest-frame UV continuum or, 
when all available data falls redward of 2 800 A, we ex- 
trapolate. These results are shown in Table [TJ column 
(6). In those cases for which only an upper limit to M* 
can be computed we also estimate a conservative lower 
limit. We do so by extrapolating a flat spectrum (/„ 
oc v°) from the reddest NIR/optical detection available 
(references in Table [TJ column 4). These lower limits are 
presented as solid bars in Figures [2j [3j and [4] 

All flux densities listed in Table [T] are corrected for 
foreground Galactic dust-extinction (see Table [JJ notes 
for the details). Conversion of the magnitudes ob- 
tained from the literature to flux densities is based on 



Fukugita et al.1 (119951) for the o ptical passba nds and on 
Tokunaea fc Vaccal (|2005h and ICohen et all (|2003f) for 



the NIR passbands. The error introduced by the as- 
sumption of these photometric systems never dominates 
over the photometric uncertainties themselves and was 
safely neglected. 

4. TOTAL STELLAR MASSES 

We infer M+ for our sample from rest-frame X-band 
flux densities. The light emitted by a galaxy in the K- 
band (e.g., the MIR photometry analyzed in this work) is 
closely related to its and thus it is a reliable estimator 
(|Glazebrook et aLl l2004). It has little sensitivity to dust 
since the majority of a galaxy's stellar population has 
moved away from the birth clouds, and because the NIR 
bands are virtually unaffected by dust extinction. Such 
derivation of M+ is more physically meaningful than an 
optical/UV luminosity; it effectively integrates over the 
accumulated M t and merger history and can only in- 
crease with time, in contrast, for instance, to UV light. 

In order to obtain M* we apply: 



M*(M ) = 2.67 x 10 _4S x 



4irD 2 L U{v ohs ) M* 



1 



L K Icst 



(1) 



where for any given object, Dl is its luminosity distance 
in cm; f v {v hs) is its flux density at the observed wave- 
length in /LtJy; observations should have been made at 
wavelengths 1.9 < v \^ s /{\+z) < 2.5 /im (e.g., this 
work); and the factor of 2.67 x 10^ 48 converts the sec- 
ond term in equation [1] to units of solar luminosity. The 
third term in equation [TJ M+/ LK Iost , also in solar units, 
must be estimated for each object. 

M+/Lif roat depends to some extent on the com- 
posit ion of the stella r population dPortinari et al.l 
20041) or, according to lLabbe et ail (|2005l ) who used 
Bruzual fc CharlotJ (|2003j) with a Salpeter IMF (jSalpeterl 
1955h . on the rest frame U — V colour, age, and 
GR B host galaxie s are blue, young , and 
[e.g. iLe Floc'h et all 120031: iBerger et all 12005 



M*. 
faint 



IChristensen et al.l 12004) . 



In ICastro Ceron et al.l (12006 

to be ~O.1M Q /L 



M*/Lic resi was assumed to be ~O.1M /L to ob- 
tain robust lower limits. For this work we compute 
M*/L Ktttat for GRBs 980703, 000210, and 000418 with 
the rest-frame if-band flux densities from Table [TJ col- 
umn (4), and M* value s derived from stellar p opula- 
tion model SED fitting (|Michalowski et all 12008) . and 
obtain the following results: O.29M Q /L for GRB 
980703, O.63M /L Q for GRB 000210, and O.43Af /L 
for GRB 000418. These results yield a mean M+/L Kiost 
= 0.45MfT)/L(7), f ully c onsistent with the average value 
in ICourtv et al.1 (12007D. and amon g the lowest ratios 
resented bv iPortinari et all (|2004D for a Salpeter IMF 
Salpeter! IT955). It is sensible to calculate an average 
M+/LK, eBt because this ratio is nearly constant, with 
little dependence on the previous SF history. In fact, 
M^/Lif rcat varies only by a factor of two between ex- 
tremel y young and extremely old galaxy stellar popula- 
tions (jGlazebrook et all 120041 ). So we estimate a con- 
servative 0AMq/Lq in the calculations of M* for our 
host sample. Table [TJ column (5) summarises our M* 
estimates. Errors quoted are statistical. We present a 
histogram of the distribution of M+ in log space f or our 
GRB host sample in Figure [TJ IVan der Wei et alj ([20061 ) 
examined redshift dependent systematics in determin- 
ing M* from broad band SEDs. They found no signif- 
icant bias for iBruzual fc Chariot! (|2003l) models with a 
Salpeter IMF (|SalDeterJll955ir ~ 

For compa ris on w e plot Figure 2 from the 
ISavaglio et alj (|2007l ) preliminary analysis in the 
background of our histogram. The two samples have 
a 25 object overlap. Our results clearly suggest more 
massive hosts, about half an order of magnitude higher 
(median M + = 10 9 - 7 M W in ou rs vs. median M+ = 
10 9 - 3 M Q in ISavaglio et al.ll2007t both distributions have 
a la dispersion of 0.8 dex, and in b oth cases the average 
has t he same value as the median) . ISavaglio et al.l (|2006l . 
120071 ) fit the optical-NIR SEDs of their host galaxy 
sample together with a complex set of SF histories. 
We can reproduce the median and average M* in 
ISavaglio et all (|2007l) with our dataset by applying 
M*/L KiBBt = 0.2 M Q /L Q (The Kolmogorov-Smirnov 
test indicates to a high probability, p ~ 0.99 likely 
becaus e of the 25 object ov erlap, that our data set and 
that of lSavaglio et al'J r2007 come from a population with 
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the same specific disctribution) . Their methodology 
therefore seems to contain an effective M^/Lk^ ratio 
around 0.2. An adjustment by a factor of ^2 to this 
ratio thus explains t he discrepancie s in M + between 
our work and that of ISavaglio et al.l (|2007[ ). We note 
that such an adjustment is within the spread of our 
calculated values (cf., 0.29 M /L , 0.43 M /L , and 
0.63M q /Lq). The fact that we find larger M* may 
also be indica t ive of underestimated dust extinction in 
ISavaglio et all (|2007t ) (see jgj}. 

Our M* are always lower than those of the normal 0.4 
< z < 2 galaxies from the Gemini Deep Deep S urvey 
fGDDS: [Abraham et al.ll200i ISavaglio et al.ll2006ri . The 
GDDS is a deep optical-NIR (K < 20.6) survey complete, 
for the already mentioned redshift range, down to M* = 
10 10 - 8 M Q for all galaxies and to M* = 10 10 - 1 M for star- 
forming galaxies. In our host sample at least 70% of the 
galaxies have M± < 10 101 M Q . This comparison clearly 
highlights the efficiency of the GRB selection technique, 
against that of traditional high-redshift surveys, to pick 
low-M* galaxies at high redshifts. 

We plot M* as a function of the redshifts for our 30 
GRB host galaxy sample in Figure O and find no intrin- 
sic correlation between the two variables. The scatter 
of M* is rather uniform across most of the redshift dis- 
tribution. Hosts with very low M+ are only found at 
low redshift. For instance, the four GRB hosts (i.e., 
060218, 060614, 030329, and 980425) with the lowest 
M* (< 10 9 M©) have some of the lowest redshift values 
in our sample. Very low-M*, high-redshift hosts would 
have been excluded since most of our largely pre-Swift 
redshifts were measured in emission, what selects pref- 
erentially bright host galaxies. Because the redshift is 
a requirement for inclusion in our sample we are effec- 
tively biased against faint systems. This situation has 
now been corrected in the Swift era when most redshifts 
are secured via afterglow absorption spectroscopy. The 
upper limits in the vertical bars of Figure [2] (i.e., the 
distribution for each non-detection measurement of M*) 
mark the sensitivity- limited curve for M*. Conversely, 
the absence of high- M*, low-redshift hosts suggests that 
such GRB host galaxies are rare. 

5. STAR-FORMATION RATES 

We compute the unobscured SFR for each host by 
means of their UV continuum luminosity. We convert 
flux densities into luminosity densities using L v (v vos t) = 
4 7 r J D|/,(^ obscrved )/(l + z) (|Hogg et al.l MM). Then 
we can calculate the unobscured SFRs by applying 
SFR(M Q yr- 1 ) = 1.4 x 1Q- 28 £ U V (erg s" 1 ) to th e rest- 
frame A = 2 800 A flux densities (|KennicuttJ 1 1998ft . The 
results are summarised in Table [U column (7). Errors 
quoted are statistic al. In addition t here are systematic 
errors of order 30% (|Kennicuttlll998h . 

The specific SFR <p = SFR/M* gives an indication of 
how intensely star- forming a galaxy is. In Figure [3] we 
plot 4> versus M* for our GRB host sample. The absence 
of hosts in the lower-left corner is explained as a combi- 
nation of selection effects and low-number statistics. A 
host galaxy in this region of the plot has both low M+ 
and low SFR, making its detection difficult unless at very 
low redshifts. As the sampled comoving volume becomes 
smaller because of the lower redshift required to make a 
detection, the chance of finding a host decreases accord- 



M K (AB) 

-16 -18 -20 -22 -24 
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Fig. 3. — Specific SFR (<f>) as a function of total stellar mass 
(M*) for our 30 GRB host galaxy sample. Squares are detections 
and triangles mark upper limits for either M t , SFR, or both. Yel- 
low diagonals: Equivalent to the vertical bars in Figure [2] Each 
yellow diagonal could be displaced vertically by the size of the cor- 
responding SFR error bar (Table [Tl column 7). Dashed diagonals: 
Constant SFRs of 100, 10, 1, 0.1, and 0.01 Moyr -1 respectively. 
Black arrows: Magnitudes of the displacements due to extinction 
by dust in the UV (e.g., 1 mag; vertical) and changes in our estima- 
tion of M* / £ K re st ~ OA Mq/ Lq (e.g., a factor of 50%; diagonal). 
Right axis: Color term equivalent to <j>. Top axis: Absolute X-band 
AB magnitude, equivalent to M*. The top/right axes represent our 
GRB host sample in color-magnitude space, effectively equivalent 
to 4> vs. M*. 

ingly. Given the size of our sample, it is reasonable to 
expect no detections in this area of the plot. The four 
GRB host galaxies with the lowest M* (GRBs 060218, 
060614, 030329, and 980425) are all at very low redshifts 
and UV bright. We also note that our sample may be 
biased against low-SFR hosts, since many redshifts have 
been measured from emission lines. On the other hand, 
the non-detection of any GRB host in the upper-right 
corner of Figure [3] should not be due to a selection ef- 
fect. Such hosts either do not exist or their afterglows 
were obscured by dust, thus preventing their localiza- 
tion. The sample at hand offers some indication as to 
the former posibility. Our two host galaxies with the 
highest M* (GRBs 030328 and 980613) would require a 
dust extinction of Ay <~ 5 mag to show there. Yet such 
dust-extinction levels can be ruled out by the constraints 
on the SFR from the IR and the radio (see $6] below). 

To exemplify how the estimation of M+/ £if reat ~ 
0.4 Mq/Lq affects the location of our hosts in the plot 
we suppose a 50% uncertainty. We repeat the exercise for 
an UV-continuum dust-extinction of 1 magnitude. The 
corresponding displacements are plotted in Figures[3]and 
[4] with black arrows. The magnitude of these displace- 
ments is limited enough not to affect our analysis. 

6. ANALYSIS AND DISCUSSION 

We find the GRB host galaxies in our current sam- 
ple possess a wide range of properties; with 10 7 M < 
M* < 10 11 M ; and 10" 2 M© yr" 1 < unobscured SFR < 
lOM yr _1 . Yet this diversity points towards low M*, 
star-forming systems. 

Part of our host sample is extingui shed by dust. 
GRB hosts 970828, 980613, a nd 990705 dLe Floc'h et"ai] 
120061 24 fim flux densities; ICastro Ceron et all 1200a 
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4. — Specific SFR (<j>) as a function of total stellar mass (M*) for our GRB host galaxy sample and other representative types 
of galaxies (for a similar plot sec Erb ct al. 2006; Castro Ccron ct al. 2006). Blue squares and triangles: Host galaxy sample from this 
work, with unobscured SFRs derived from the rest-frame UV continuum. Triangles mark upper limits for either M*, SFR, or both. The 
open blue triangle and squ are mark the hosts of GR Bs 970828 and 980613 respectively (see Black squares: SFR values constrained 
with a best fit SED model (Castro Ccron ct al. 2006). Grey squares: Highest/lowest SFR values for those hosts for which a best fit model 
could not be established (Castro Ccron ct al. 2006). Both the black squares and the grey squares have been shifted here corresponding ly 
to compensate for the difference in M 4 /Lj{ rB!t methodology (i.e., from a lower limit value of 0.1 Mq/ Lq in Castro Ccron ct al. 2006 to 
a best estimated value of 0.4 Mq / Lq in this work). Yellow diagonals, dashed diagonals, and black arrows are as in Figure[3] Right axis: 
SFR timescale (TgpR = M*/SFR), the inverse of </>. On this scale the solid horizontal lines represent the age of the universe for the marked 
redshift. The distribution of our sample i n parameter space suggests that GRBs tr ace galaxies that are not se lected with oth er techniques. 
Data p oints: GRBs f rom t his w ork and ICastro Ceron et al.l l|2006l) . DRGs from Ivan Dokkum et all (120041') . LAEs from [Gawiser et al.l 
pOO??). INilsson et all S200% . and ILai et al.l (|2008I 1 (filled stars: average values of the LAE populatio n as a whole, obtain ed from stacked 
photometry); N. Pirzkal, S. Malhotra, J. E. Rhoads, & C. Xu (2007, priv. comm. ; empty stars; see | Pirzkal et al ] l|2007| ) for the average 
value s and a desc r iption of these s ources). Spectroscopica lly confirmed LBGs from Shapley et al. (2001) and Barmby ct al. (2004). SMGs 
from I Borvs et aLl 120051 . M+) and ICha pman et al.l ((2005 , SFR), where we have considered Lbol — ^farlR; then applied the Kcnnicutt 
(119981) calibration, z ~ 2 from lErb et al.l ( 120031" ) and lReddv et al.l l Hffl06T) . 



SEP fitting), as well as 980703, 000210, and 000418 
(ICastro Ceron et al.ll2(M SED fitting for GRB 980703; 
iMichalowski et all 120081 SED modelling), have highly 
obscured SFRs. Additionally several authors argued 
for dust extinction in the host of GRB 31203 (e.g., 
iProchaska et al J 12004: iMargutti eFail 120071) . Applying 
the recipe in ICastro Ceron et aLl (|2006l ) to this host 
galaxy's MIR photometry (/3.6 pm = 216 /iJy ± 3 /xJy; 
f5.8fj.rn = 390 ^Jy ± 16/xJy; /24 pm = 13103/xJy ± 
41 ^Jy; flux de nsities corrected for foregroun d Galactic 
dust-extinction. lLutz| [r999: Dutra ct al. 2003) we obtain 
a SFRl 8 _ 1000 = 13 Mq yr _1 . That brings the total num- 
ber of extinguished hosts to at least 7 out of 30, and 



allows us to crudely estimate that >25% of the sam- 
ple in this work suffer significant dust extinction (Ay > 
lmag). Neither our host galaxy sample, nor those oth- 
ers cited in this work, are bias-free. The searches for 
the GRBs in such samples have been carried out mostly 
following the localization of an optical afterglow, implic- 
itly biasing the sample against dust-extincted systems. 
Such potential bias strengthens our statement on dust ex- 
tinction in GRB host galaxies. Parenthetically, we note 
that GRBs 970828, 980613, 980703, and 990705 make 
up two-thirds of a redshift-z^l-selected, small sub sample 
(|Castro Ceron et al.|[200l iLe Flocm et al.ll2006f) . They 
hint at the possibility that even a higher fraction of hosts 
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are affected by dust extinction, though with the caveat 
o f low-number statist i cs. 

ICastro Ceron et all (|2006l ) plotted <j> versus for six 
GRB hosts and samples of five other representative types 
of galaxies: distant red galaxies (DRGs), Lya emitters 
(LAEs), Lyman break galaxies (LBGs), submillimeter 
galaxies (SMGs) and an ensemble of optically-selected, z 
~ 2 galaxies from the Great Observatories Origins Deep 
Survey-North field. In Figure [4] we plot, with blue sym- 
bols, our 30 host sample using a revised M^/Lk^ r atio, 
along with the samples in lCastro Ceron et alj (|2006f ) ■ In 
our (f> versus M* plot (i.e., Figured) we use SFRuv) the 
obscuration of SF by dust pulls the GRB data points 
down along a vertical li ne. One way to reconcile t he <f> 
values of host galaxies in ICastro Ceron et"aTI (|2006f) and 
this work is to invoke extinction by dust of the order of 
Ay ~ 1-3 mag (see below ) . The conversion from A\jy to 
A v follows ICardelli et all (fl989h . 

A primary scientific goal in the quantification of galac- 
tic evolution is the derivation of the SF histories, as de- 
scribed by th e temporal evolution of th e star-formation 
rate, SFR(i). ICastro Ceron et all (|2006l ) noted that their 
sample had Tsfr < iunivcrsc, allowing for a history of con- 
stant SF, with a robust lower limit in M* (M*/ L,K ICBt ~ 
0.1Mq/L q ). For the sample we present in this work, 
where we adopt M+/L}( T<xrt ~ 0.4 Mq/Lq, clearly few 
GRB host galaxies are not allowed to have a history of 
constant SF (i.e., young stars dominating the stellar pop- 
ulations of old galaxies; see the right ordinate axis in Fig- 
ured]). Either a starburst episode was present in the past, 
or a higher recent SFR is required. The latter possibility 
is consistent with a fraction of GRB hosts having SF ob- 
scured by dust. The hosts of GRBs 970828 and 980613 
(open blue symbols in Figure |4|) are good examples be- 
cause, under the assumption of constant SF, major dust- 
extinction must be invoked to account for the age differ- 
ences. 0uv estimates result in Tsfr ~ 6 Gyr for GRB 
970828 and Ts fr ~ 32 Gyr for GRB 980 613, while cj> Llc 
estimates (see ICastro Ceron et all I2006D result in Tsfr 
~ 300 Myr for both of them. The discrepancies in Tsfr 
imply a dust extinction of the order of Ay ~ 1.6 mag for 
GRB 970828 and A v ~ 2.5 mag for GRB 980613. These 
discrepancies are consistent with the radio-constrained 
SFR upper limits (-TOOMgyr" 1 for GRB 970828, and 
~500M^yr-^ for GRB 980613) derived by applying the 
lYun fc Carillil (|2002D me thodology to the deepest radio 
upper limits reported bv lFrail et aU (|2003h . 

A dilution effect is present in our MIR photometry. 
Hosts in our current sample are not spatially resolved in 
the Spitzer imagery (in the case of GRB 980425 we utilize 
the total flux of the galaxy for consistency with the rest 
of the sample). To estimate their M* we measure the 
total if-band light. Lk traces the accumulation of M* 
([Glazebrook et aLll2004f ) while, most commonly, the SF 
is ongoing in only a small part of the host galaxy. So we 
do not normalize our sample's unobscured SFRs by the 
total stellar mass of the star-forming region(s), rather 
by M+, which results in lower <f> values. This dilution 
effect pulls the GRB data points in a versus M* plot 
down along the diagonal (dashed) lines marking constant 
SFRs. 

Extinction by dust, coupled with the dilution effect, 
could be responsible for the apparent envelope that 
can be visualized in Figures [3] and [4] a flat plateau 



with no objects above a certain (f> value (~2.5 Gyr -1 ), 
and that starts to curve down beyond a particular M* 
(10 10 Af Q ). Correcting for dilution and, chiefly, for dust 
extinction would yield a new plot where our host sample 
would align consistently w ith the results/upper limits of 
ICastro Ceron et "all (|2006| ). and provide support to the 
claim that GRB host galaxies are small and have some 
of the highest <j> values. 

We conclude by putting forward a simple idea for GRB 
hosts based on the data analyzed here. As a working 
hypothesis we suggest that, while low M* hosts might 
contain no dust (i.e., host galaxies with a low M+ and a 
low SFR are rare; see $5} , progressing upwards in the M+ 
distribution of host galaxies will yield significant dust ex- 
tinction, as well as the already mentioned dilution effect 
(i.e., the apparent envelope described above for Figures[3] 
and Our suggestion is consisten t with the theoretical 
predictions presented in lLapi et a l. (2008). They further 
predict that GRB host galaxies trace the faint end of the 
luminosity function of LBGs and LAEs. Future work 
(i.e., Herschel observations) on a complete host sample 
will allow us to test this by quantifying dust extinction 
and the dilution effect. 

The nature of GRBs 060505 and 060614 is strongly 
debated as no supernova was associated with these 
long- d uration GRBs to d eep limits (iFvnbo et alJ 



20061: lOehrels et all 120061 iDella Valle et all 12006a! : 
Gal- Yam et al.l I2006D . GRB 060505 falls within the 
distribution of other long-duration GRB hosts in our 
sample, whereas GRB 060614 seems to be an outlier. 
Though this may indirectly suggest that the progenitor of 
GRB 060614 is different from other typical long-duration 
GRBs, we note that its SFR is in range with that of the 
bulk of the sample; and as for M*, its properties are not 
very different from some of our other low-redshift host 
galaxies (e.g., GRBs 060218, 030329, and 980425). 
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